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Abstract 

An effect of polarization-mode dispersion, nonlinearity and random variation of dispersion 
along an optical fiber on a pulse propagation in a randomly birefringent dispersion-shifted optical 
fiber with zero average dispersion is studied. An averaged pulse width is shown analytically to 
diffuse with propagation distance for arbitrary strong pulse amplitude. It is found that optical 
fiber nonlinearity can not change qualitatively a diffusion of pulse width but can only modify a 
diffusion law which means that a root mean square pulse width grows at least as a linear function 
of the propagation distance. 
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Polarization-mode dispersion (PMD), which is a pulse broadening caused by random 
variation of optical fiber birefringence, has recently become a major drawback in the devel- 
opment of new high-bit-rate optical communication systems Another 
effect, which limits bit-rate capacity, is pulse broadening caused by group- velocity dispersion 
(GVD). Use of a dispersion-shifted fiber with zero average GVD can reduce this effect, how- 
ever, in such fibers GVD inevitably fluctuates around zero along the propagation direction 



and hence pulse broadening still occurs 
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Nonlinearity in optical fibers results in 



the coupling of both PMD and GVD effects, so in general they can not be studied separately 
in contrast to linear case. Linear PMD was first studied in Refs. while nonlinear 



PMD was addressed in numerical experiments |2j and analytical studies based on a pertur- 



bation expansions around soliton solutions of deterministic equations \c , 
effect of random variation of GVD was studied in Refs 
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Here an exact 



analytical (nonperturbartive) theory is developed for the case of fiber with random birefrin- 
gence and random GVD with zero mean and arbitrarily strong nonlinearity (arbitrary pulse 
amplitude). No assumption like closeness to any type of soliton solution is necessary for the 
results of this Article to be valid. The main result is that a statistical average (over random 
variation of fiber parameters) of root mean square pulse width Trms grows with distance 
at least as a linear function of propagation distance. This means that random diffusion of 
optical pulse width can not be prevented by an arbitrarily strong nonlinearity. It is shown 
that random diffusion fundamentally limits the bit-rate capacity of an optical fiber. 

Neglecting second-order GVD (dispersion slope) effects, stimulated Raman scattering 
and Brillouin scattering, the propagation of optical pulses in birefringent optical fibers is 
described by the two-component vector nonlinear Schrodinger equation (VNLS) p. Il4 . Il9 . 2(] ] 
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idzVa + ^2 ^c(3{z)rip + i ^ rn a(3 {z)d t rj P + d(z)dfr] a 



(3=1 



0=1 



+a{z)N a {r 1 ) = iG{z)r] 



(1) 



where z is the propagation distance along an optical fiber, 771 and r] 2 correspond to the 
complex amplitudes of two orthogonal linear polarizations, t = r — zjc\ is the retarded time 
and r is the physical time, q is the speed of light, and d(z) is the dispersion, which is related 
to first-order GVD /3 2 as d(z) = —\(ih,(z). The right hand side (rhs) of Eq. (JTJ) describes 
linear losses and amplifiers, G(z) = ^— 7+[exp (z a -y) — l\T,^ =1 5(z— z k ) \ , a = (27m 2 ) / (XoA e f f) 
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is the nonlinear coefficient, n 2 is the nonlinear refractive index, Ao = 1.55/im is the carrier 
wavelength, A e ff is the effective fiber = kz a (k = 1, . . . , N) are the amplifier 

locations, z a is the amplifier spacing, and 7 is the loss coefficient. Distributed amplification 
can be also included by adding ^-dependence into 7. Properties of fiber can be different 
along optical line, e.g. A e ff could be different if line consists of several pieces of fiber with 
different cross section, and, respectively, coefficient a generally depends on z. In a similar 
way, all parameters of fiber, like d(z) also depend on z. 

The self-conjugated matrices A(z) and fh(z) describe, respectively, the differences in wave 
vectors and the anisotropy of the group velocities of the two modes corresponding to the two 
different polarizations. Both matrces A and fh are made traceless. The trace of the matrix 
A is excluded by a phase transformation f) — > 77 exp(i<poz) . The trace of the matrix fh is 
zero because Eq. Q is written in a frame moving with average group velocity (note that 
group velocity is generally ^-dependent). It is assumed in Eq. that the dispersion d(z) 
and nonlinearity are isotropic because their anisotropy is usually negligible in optical fibers. 
Vector N = (Ni, N2) T , which represents the contribution of Kerr nonlinearity, is given by: 

= [(|^| 2 + ^ 2 | 2 )^l + , 
iV 2 (*) = [(l\*l\ 2 + \*2\ 2 )*2 + 1*1**2] (2) 

(see 0,0). 

The change of variables £ = rje" ^0 G i z ') dz ' ( see e .g. Refs. |2ll l22|) removes rhs of Eq. (JT|) 
and gives: 
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+d(z)d^ a + c(z)N a (O = 0, (3) 

where c(z) = cr(2;)exp ^2 J Q 2 G(z')dz'\ Thus all linear fiber losses and amplifications are 
included into coefficient c(z). 

The isotropic case, which corresponds to zero matrices A = fh = 6, allows a solution of 
Eq. © with constant polarization, e.g. £1 7^ 0,^2 = 0. Components of matrices A and fh 
fluctuate strongly as functions of distance z. Fluctuations correspond to violation of circular 
symmetry of the fiber. The matrices A and fh change in optical fibers with time on a scale of 
few hours because of environmental fluctuations, however, for typical optical pulse duration 



(lOps), one can consider A and m as functions of z only. It means that disorder is frozen in 
the fiber. The matrix A gives the leading order contribution in Eq. © because a typical 
beat length Zb ea t (typical length at which a relative phase shift between two polarizations 



12 1 . Contribution of all other terms 



caused by A becomes ~ ir) is ~ 20m in optical fibers 
in Eq. Q are essential on a scale of the order of 10km and larger. Thus, it is convenient to 
introduce a slow variable \l/ as 

2 

(3=1 

where the unitary matrix T(z) is given by the solution of the matrix equation: 

2 

:^A Q<5 (z)T 5/3 , 2^(0) = J, (5) 



dz 

5=1 

where I is the identity matrix. 

The slow variable ^ changes on the scale ~ 10km because all fast dependence on A is 
included into t. Thus E q . ffl ca_n beaded over distance caueh ,ar g e r than , bm „ but 



still much smaller than 10km 3|. 



averaged equation 



1 HQ: 



Here we use the simplest possible form of an 
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0*tf « + i m a p{z)d t ^p + d{z)d^ a + N a = 0, (6) 
P=i 

where m(z) = T _1 ■ m(z) ■ T and the nonlinear terms are 

N 1 = c(z)(\y 1 \ 2 + ^ 2 \ 2 )V U 

N 2 = c(z)(^ 1 \ 2 + \y 2 \ 2 )V 2 . (7) 

The typical scale z corr of variation of m in optical fibers is between 10m and 100m. Typical 
scales at which a pulse experiences essential distortion are a dispersion length z disp = t^/di, 
d\ is the typical amplitude of dispersion variations, to is the typical pulse width; a nonlinear 
length z ni = l/(c(z)p 2 ), p is the typical pulse amplitude; and PMD length z m = t /m , 
mo is the typical amplitude of variation of matrix m components. As a typical example 
one can set to ~ lOps, d\ ~ lps 2 /km, p 2 ~ 2mW,c(z) ~ 0.001 (km mW) -1 , and mo ~ 
lps/km. One gets Zdi sp ~ 100km, z n i ~ 500km, and z m ~ 10km. Thus, a minimal length 
of pulse distortion z pu i se is ~ z m ~ lOfcm. the length z pu i se is much larger than z corr and, 



according to the central limit theorem, rh can be treated at scales of the order of z vu \ se as 
random Gaussian processes with zero correlation length and zero mean (rh) = 6 ((. . .) means 
ensemble averaging over statistics of rh(z) and d(z)). 

The traceless matrix rh can be represented in terms of Pauli matrices: rh(z) = 
Sj=i m j( z )&j> where the correlation functions of components of the real vector m are given 
by 

(m j (z 1 )m k (z 2 )) = M j 5 jk S(z 1 - z 2 ). (8) 

The vector M does not depend on z and M is defined from the original problem with short 
but nonzero correlation lengths as Mj = J (mj(z)mj(z'))dz' . 

Similarly, dispersion d(z) in a dispersion-shifted fiber (where average GVD is shifted to 
zero but random variations of GVD are essential Q, S, Q) 

can be represented by random 
Gaussian process with zero correlation length and zero mean: 



(d{z 1 )d{z 2 )) = D5{ Zl - z 2 ), (d{z)) = 0, (9) 

where D — J (d(z)d(z'))dz f . The quantities rh(z) and d(z) are assumed to be the independent 
random processes: (mj(z)d(z')) = 0. 

Below the Eqs. © — are used to obtain exact result on evolution of pulse width along 
z. First one can conclude that as (fh(z)) = and (d(z)) = 0, there is no preferred direction 
along t and the average pulse position is zero: 

(t) = (Jt(\y 1 \ 2 + \V\ 2 )dt)/P = 0, (10) 

where 

P = J(\^ 1 \ 2 + \^ 2 \ 2 )dt (11) 

is the time-averaged optical power. The quantity P is an integral of motion of VNLS (JHJ) : 
d z P z = 0. 

Consider now evolution along z of 

A = j(t-(t)) 2 (\^ l \ 2 + \^\ 2 )dt, (12) 

which is related to the root mean square pulse width T RMS by 

T 2 RMS = A/P (13) 



Using Eqs. (0) and (fTUJl . integrating by parts over t, and applying vanishing boundary 
conditions at infinity one derives 
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A z = d{z)B {d) -J2 m A z ) B{mj \ ( u ) 

where 

J a=l 

B {mi) = -2 J t(¥jtt 2 + c.c.)dt, 

B {m2) = 2% J t{^ 2 ~ c.c.)dt, (15) 
B {m3) = -2 J t(|^i| 2 - \^ 2 \ 2 )dt. 

It is essential that all expressions in rhs of (|T3|) do not have explicit dependence on random 
variables d and m, which allows one to differentiate them over z and find, using again Eq. 
© and integrating by parts, that 

#W = 8d(z)X -2c(z)Y + 0(m), (16a) 
B (mi) = _ 2mi P + 0(d, m 2 , m 3 , c(z)), (16b) 

£(ma) = _ 2m2 p + mi, 77l 3 , c(*)), (16c) 

S (ms) = -2m 3 P + 0(d,m h m 2 ), (16d) 



where 



F = / (|^| 4 + |^ 2 | 4 + ^xH^I 2 )^. (17) 
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Notation 0(d,m2,m 3 ,c(z)) in Eq. (|16b|) means extra terms which are linear in 
d,m2,rri3,c(z). And the same type of notation is used for the similar terms in rhs of 
Eqs. (|16a|) . (|16c|) . (|16d|) . Explicit expressions for these terms are not given here because 
they are bulky and, as it is shown below, they vanish after statistical averaging. Note that 
or m = 0, ^2 = and d = Const Eq. (j!6a)l coincides with a so-called virial theorem 



3, 



24, 



25| . However, direct application of the virial theorem to Eq. (JBJ) is not possible 



because it would require determination of A zz by differentiating Eq. (JT3j) over z, which 



would result in appearance of terms including derivatives of random variables d(z) and m 
over z. Here these problems are avoided by studying B z d ^ and B z m ^ instead of A zz , but 
as a result of such procedure, one obtains cumbersome expressions (|16a|) — (jl6d|) compare 
with t 
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24 



re compact expression for the virial theorem in deterministic scalar case (see Refs. 

3)- 



Consider the statistical average of 

= d{z) (b® (0) + jT Bf (z')dz^ , (18) 

and substitute B^ (z') in that equation by rhs of Eq. (jlfiaj) . Assume for a moment (to choose 
a correct limit z corr — > 0) that the correlation length z corr of all random processes d(z) and 
m(z) is small but finite and that all correlation functions decay at least exponentially or 
faster with distance. Now choose a distance zq in such a way that z corr <C zq <C z pu i se . The 
statistical average of Eq. (fTS|). using Eq. ()16aj) . gives 

(4 d) ) = (d(z) f (u{z')X{z') - 2c(z')Y(z'))dz'), (19) 

where all corrections to that equation are exponentially small and vanish in the limit z corr — > 
due to casuality constrain. Casuality constrain implies here that ^(z) does not depend on 
d(z) and m(z) if z < z. 

The condition z <C z pu i se allows one to write solution of VNLS © as ^(z') = ^(zo) + 
0(z' — Zq) for z < z' < z. Taking the limit z corr — *> 0, for which z — > z, and using Eq. ()19j) 
one obtains 

(A { z d) ) = AD(X). (20) 
To average Eq. (j!4)l one rewrites it in the equivalent form 

A z = d{z) (BW>(0) + jT B { f(z')dz'^ 

-J2^A z )( Bimj) (°)+ f B { p\z')dz'^, (21) 

and substitute expressions for B z d \ B z mi \ Bi m2 \ Bz from rhs of system ()16a|) — (|16d|) into 
Eq. (|21j) . As a result, one finds, by averaging over distribution of m and d, that 

3 

<A,> = PX>i +4D <*>- (22) 



Here (P) = P because P is an integral o motion. All terms in rhs of Eq. (|22|) are positive 
definite and hence (A) grows with z. 

Using Eqs. f|13jl and ()22|) one obtains an expression for the statistical average of the root 
mean square pulse width 

3 AD f z 

(Trms(z) 2 ) = T RMS (0) 2 + zVM i + - (X(z'))dz'. (23) 

3=1 J ° 

Eq. (|2*Hjl is exact for system (jHJ) — (JSJ) because Eq. (|2*Hjl is derived for arbitrary strong 
nonlinearity. Thus, this is essentially nonperturbative result. Remarkably, Eq. ()23|) does 
not explicitly depend on the nonlinear coefficient c(z). This fact is very peculiar propertiy of 
VNLS © and is related both to the generalization of the virial theorem of scalar nonlinear 



Schrodinger equation [2j 



24. 



and to the zero correlation length limit of fiber parameters 
(JHJ) and (JHJ). In the linear case, which corresponds to c(z) = 0, X does not depend on z 
and the growth of (T R ms( z ) 2 ) with z is linear (i.e. the pulse width experiences diffusive 
growth caused by random variations of GVD and PMD matrix rh(z)). For nonzero c(z), 
the nonlinearity is only responsible for nontrivial dependence of (X) on z, and therefore, 
for modification of the (still) diffusive law. The term P X^=i Mj is constant for arbitrarily 
strong nonlinearity which means that growth of the pulse width can not be slower than linear 
in distance z. This gives a fundamental limit to the bit-rate of information transmission in 
optical fiber systems. 

Estimating X, P, |M| and D by p 2 /t , p 2 t , z corr ml and z corr d\, respectively, one derives 
from Eq. © that (3 = ((T RMS (z) 2 ) - T RMS (0) 2 ) /T RMS (0) 2 

~ zz corr (m 2 ) ft 2 , + d\/t^\. Then the minimal requirement for small information loss, (3 ~ 1, 
results in a limitation for the pulse width: to ~ lOps, for z ~ 10 3 A;m, and the typical values 
Zcorr ~ 100m, mo ~ lps/km, d\ ~ lps 2 /km. It suggests that construction of high-bit-rate 
lines based on the dispersion-shifted fiber requires essential improvement of fiber production 
technology, and/or implementation of both PMD compensation j3, [22 1 and the pinning 
method [lpj, to reduce both PMD and random GVD variation effects. 

Note that the extension of the results of this Article to the case of both nonzero average 
dispersion (dispersion-shifted fiber with nonzero average dispersion) and dispersion-managed 
systems is an open problem. E.g. nonzero (d) = do(z) would result in appearance of new 
terms 8d (z) f z [d {z') {X (z')) - 2c(z'){Y(z'))]dz' in rhs of Eq. (j22) which generally are not 
sign-definite. In the case of dispersion management these terms oscillate fast with distance 
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z so their contribution to Eq. (|22|) can be small. Additional research is necessary for the 
case of nonzero average dispersion which is outside the scope of this Article. 

In conclusion, the exact analytical expression (|23|) for random diffusion of the averaged 
optical pulse width is derived. Eq. (|23|) introduces a fundamental limit on the minimal pulse 
width for which information transmission is possible in nonlinear dispersion-shifted optical 
fiber systems with zero average dispersion. 
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